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Abstract 

We discuss magnetic Schrodinger operators perturbed by measures from the 
generalized Kato class. Using an explicit Krein-like formula for their resolvent, 
we prove that these operators can be approximated in the strong resolvent sense 
by magnetic Schrodinger operators with point potentials. Since the spectral 
problem of the latter operators is solvable, one in fact gets an alternative way 
to calculate discrete spectra; we illustrate it by numerical calculations in the 
case when the potential is supported by a circle. 

1 Introduction 

Schrodinger operators are used for modelling a particle confined in a quantum- 
mechanical system. Depending on a potential which describes how the particle inter- 
acts with its environment, one can consider a wide range of physical situations. In 
this paper, we are particularly interested in potentials in dimension two supported 
by zero measure sets; the supports could be for example graphs, curves or points. 

The motivation to study such operators is based on the fact that they represent 
simple mathematical models of various nano-structures like quantum wires, photonic 
crystals, quantum dots, etc. One possible way to describe them is via quantum 
graphs; it means that one considers ordinary differential equations on the graph 
edges, which are coupled through boundary conditions at the graph vertices so that 
the resulting operator is self-adjoint, see jKj or |KSj . The operators we are going to 
deal with yield an alternative approach. The particle is not confined to the graph, 
but it moves in its vicinity if the potentials are attractive. Hence the latter model is 
in a sense more realistic and it enables us to take the tunnelling effect into account. 

We aim to prove a limit relation between two classes of operators in L 2 (M 2 ) in the 
presence of a magnetic field: those with attractive potentials supported by a curve or 
a graph on one side, and operators with point potentials on the other side. The crucial 
feature of the latter operators is the solvability of their spectral problem; if the number 
of potentials is finite, then the essential spectrum stays unchanged and the discrete 
spectrum can be calculated numerically by solving an implicit equation. Therefore, 



if we are able to find a sequence of point potential operators which approximates the 
given operator, we get an approximate method to calculate its discrete spectrum. 

In fact, we will show that the approximation works for a larger family of po- 
tentials than just the ones supported by a curve. The regular potentials with Kato 
property will be also included; that is why we speak about a generalized Kato class 
and potentials are replaced by more general measures. For a given measure m from 
the generalized Kato class, it is possible to define the operator —A + m via asso- 
ciation with a closed and semi-bounded quadratic form, see e.g. jBEKS and |SV| . 



The second way to define Schrodinger operators with potentials supported by zero 
measure sets is by prescribing the operator domains, see |Polj for general singular 
perturbations, and |AGH Hj or GHS] for point potentials; there the domains are given 
by imposing a boundary condition on wavef unctions. 

It was shown in |BFTj that the free Laplacian perturbed by a measure with 
Kato property can be approximated by point potential operators; in dimension one 
the convergence is in the norm resolvent sense, while in dimension three it is in the 
strong resolvent sense. According to |EN2j , the situation in dimension two is similar to 
the three-dimensional case, moreover, the authors presented several physical systems 
where the approximation is useful in spectral calculations. 

The present task is to prove that the approximation also works in the presence of 
a magnetic field. It turns out that the main difficulty is not the proof itself (it easily 
carries over from the non-magnetic case) but rather the lack of information about 
magnetic systems. Namely, we first need to clarify the definition of perturbations by 
a measure in sectional Then in section 0] we der ive an e xplicit formula for resolvents; 
that must be done without using results from BEKS directly because their proof 
relies on the positivity preserving property of the free Laplacian. Section ideals with 
point potentials in the presence of a magnetic field. Finally, in section El we state the 
main approximation claim and we apply the approximation to a simple example in 
section [3 where the magnetic field is homogeneous and the potential is supported a 
circle. 



2 Magnetic Schrodinger operator in R 2 

The free magnetic Schrodinger operator on L 2 (IR 2 ) is given by 

(-iV-A(x)) 2 , 

where A(x) is a vector potential, whose components A\ and A 2 belong to C 00 ^ 2 ). 
According to |CFKS| chapter 1.3], there exists a closed and positive quadratic form 
h, 

D(h) = {tp G L 2 (R 2 ) : (dj - iAj)ip G L 2 {R 2 ), j = 1,2} 

2 
3=1 
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We define the free magnetic Schrodinger operator H as the unique self-adjoint op- 
erator associated with the form h, i.e. 



D{H ) C D(h) 



ip,il>eD{H ). 



Moreover, from [CFKS[ theorem 1.13] we know that C^°(M 2 ) is a form core of Hq. 

By [BGP1 , theorem 14] the resolvent {Hq — z)~ l has an integral kernel Gq(x, y; z) 
which is continuous away from the diagonal x = y. The singularity of Go(x,y;z) 
on the diagonal is of the same type as the one for the non-magnetic Green function; 
following [BGP2] theorem 15], it can be rewritten as 



G (x,y;z) = - 

where G r en (-;z) is continuous on 
regularized Green function 



-— H\x-y\) + G™(x,y;z), 
Lix 



I 2 x R 2 . 



Hence it is possible to introduce the 



f(a; z) 



lim 

\x— al^O 



G (x,a; z) + — ln(|x - y\ 



G r en (a,a; z). 



(1) 



We will need this function when defining a perturbation by point potentials. 

In the special case of a homogeneous magnetic field B one can write Green function 
Gq(z) explicitly. For example, in the symmetric gauge, A(x) = (— \Bxz, \Bx\), 
B G JR., Green function has the following form, see |DMMj 



G (x,y;z) 



4tt 



$B(x,y)T 



B\ 



2\B\ 



U 



\B\ 



2\B\ 



\B\ 



x-y\ 



(2) 



U is the irregular confluent hypergeometric function 
factor 



13.1.33] and $b is a phase 



$B{x,y) = exp 



LB 



\B\ 



— {xiy 2 -x 2 y 1 ) ~-£-\x 



3 Perturbation by a measure 

Next, we perturb the magnetic Hamiltonian Ho by a measure —7m in the following 

way, 

H im = H - 7m, 

where m is a finite positive measure from generalized Kato class, which means in 
dimension two that it satisfies 

lim sup / I ln(|x — y |)| m(dy) = 0, 
£ ^°xm 2 Jb(x,e) 
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with B(x,e) denoting the circle of radius e centred at x. 7 is a bounded and con- 
tinuous function mapping F := supp(m) into R + , thus we consider only attractive 
potentials. An example of such measure is the Dirac measure, supported by a curve 
or graph; one can easily check that the condition above holds. 

By [SV. theorem 3.1] the potential generated by m is (— A)-form bounded with 
infinitesimally small relative bound. In order to define H im properly, we need a 
similar form-boundedness with respect to H$. We cannot use the mentioned result 
directly as it was formulated only for Dirichlet quadratic forms, i.e. the ones which 
are positivity preserving. Instead, we can employ the diamagnetic inequality to pass 
from the non-magnetic system to the magnetic one. In the following, || • || P)9 denotes 
the norm of an operator acting from L P (R 2 ) to L q (M. 2 ), 1 < p, q < 00. 

Lemma 1. Let m be a positive measure from the generalized Kato class w.r.t. —A 
and let H Q be the self-adjoint operator defined above. Then for each a > there exists 
b G K such that the following inequality holds for any i/j G Cg°(K 2 ) 

/ \iP(x)\ 2 m(dx) < ah(ip,ip) + b\\ip\\ 2 2 . (3) 
Jr 2 

Proof: Since m belongs to the generalized Kato class w.r.t. —A, by jS V|. theorem 
3.1], the following inequality is fulfilled for any ip G C^°(1R 2 ) 

/ |^(x)| 2 m(rfa;)<^[||V^||^ + Cll^ll2], (4) 
Jr 2 

where according to |SV| remark 1.7(b)], u := ||(— A + C)~ lm l|oo and thus it decays 
with ( growing to +00. It was proved by approximating m by a sequence of non- 
negative potentials V n G L 2 (R 2 ) n L°°(M 2 ) such that 

lim / \i)(x)\ 2 V n (x)dx = [ \i){x)\ 2 m{dx) V^gC °°(M 2 ). 
n ~ fCO Jr 2 Jr 2 

|SV| theorem 2.1] states that such sequence exists, ||(— A + C) _1 Ki||oo,oo < uj, n E N 
and inequality (pQ) holds also when m(dx) is replaced by V n (x) dx. 
Let us write 

/ \^x)\ 2 V n (x)dx = \\V*(H + C)-HHo + Oh\\l 
Jr 2 

< iik 1 (H + cr* \\i, 2 [hu>, w + ccv, n • 

Repeating the proof of AHS, theorem 2.5], we make use of the diamagnetic inequality 
[HSUJ in this form 

\e-tH ip\ < e-t(-A)\if)\ t > 0, V e L 2 (R 2 ). 
Then from the expression 

(Ho + C)"^ = j^iy J rh-tCe-tHo dt 



4 



we get \(Hq + C) 2-0| < (—A + 2 \ip\ } which in turn yields 

(#o + cr^ii2, 2 < ii^?(-A + o-'iii,2< n^i(-A + c)-vin 2i2 . 

The Stein interpolation theorem |RS[ theorem IX.21] and the duality between || • j^i 
and || • ||oo i00 imply 

\\Vn 2 (H + O'^Wh < IIK(-A + O'ltlK-A + C) _1 K||ioc < u. 

Finally, the convergence of V n to m and the fact that u — > as ( — > 00 finish the 
proof. □ 
Since C^°(IR 2 ) is dense in D(h) it is possible to define linear operator I m 

I m : D(h) (-> L 2 (m) := L 2 (R 2 ,m) 



I m ip = ip Vip e C { 



o°°(K 2 ) 



Then (J3J) can be extended to whole D(h) with function ip on the lhs being replaced 
by I m ip and thus I m is bounded. Now, consider quadratic form h im given by 

D(h 7m ) = D(h) 

h jm (ip,tp) = / (V$(x)+iA(x)$(x)).(V<p(x)-iA(x)<p(x))dx 



I m ip(x) I m ip(x) j(x) m(dx), 

We employ the KLMN theorem, see fRSj theorem X.17], to conclude that h im is 
lower semi-bounded and closed. Thus there exists a unique self-adjoint operator H jm 
associated with this form. 

The definition that we have presented above applies to both regular potentials 
m(dx) = V(x)dx and potentials supported by zero measure sets T. In the latter case 
there is an alternative way to define operator H im via boundary conditions. Consider 
an operator which behaves as H away from the set T 

H im ip(x) = (— iV - A) 2 ip(x) xeM. 2 \T 

with the domain consisting of functions ip such that H ip e L 2 (IR 2 \r), their restriction 
to R 2 \ T is smooth and which are moreover continuous at T and have a jump in the 
normal (w.r.t. curve T) derivatives, 

dip dip 

(x) — — — (2) = —7(2) ip{x) , x G r . 



dn + dn_ 

One can check that H im is e.s.a. and by Green's formula we have (H 1Tn ip,g) = 
h jm (ip,g) for all ip e D(H im ) and g e C °°(M 2 ). Since C^°(1R ) is a core of /i 7m , the 
closure of H im can be identified with H im . This definition is applicable to curves 
T which do not have any cusps and only a finite number of smooth edges meet in a 
node. 
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4 Krein-like formula 



Throughout this work, a crucial role will be played by Krein's formula which gives 
us an explicit expression for the resolvent of H im . Originally, the formula was used 
for Hamiltonians perturbed by a finite number of point interactions, later, it was 
generalized t o a larg e family of operators, see e.g. |Polj . 

In paper BEKS the authors derived the resolvent for the free Laplacian perturbed 
by measure m from generalized Kato class, using the positivity property of Green 
function of the Laplacian. Although their proof does not apply to magnetic systems, 
one would still expect that the resolvent corresponding to the operator H with 
perturbed by m should look the same in the presence of a magnetic field: 



R{z) = Ro(z) + R d x,m{z) 



1 

Rfn.m \Z) 
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Rm,dx(z), (5) 



where R^ u (z) is an integral operator acting from L 2 (M?, v) to L 2 (M 2 ,/i), fi and v are 
two arbitrary positive Radon measures and 

Rn, v {z)il){x) = J G (x, y; z) ip(y) v{dy) fi - a.e. 

Note that R m ,dx{z) = I m Ro(z)- 

To prove that (0) is indeed the resolvent of H^ m , we first show several auxiliary 
results. 

Lemma 2. Assume that measure m is finite, i.e. J m(dx) = l m < oo and z G p(H ). 
Then R d x, m (z) = (R m ,dx(z))* . 

Proof: For z G p(H ) the operator R m ,dx{.z) is bounded and so is its adjoint. We 
have to check that in the expression 



I = (f,R m ,d x (z)ip) L 2 {m) = J m(dy) J dx f(y) G (y,x; z)ip(x) 

one can interchange the order of integration. This is possible if the following integral 
is finite 



I<h'= J J m(dy)dx\f(y)\\G Q (y,x;z)\\'il;(x)\. 

By Kato's inequality (sec [BGP1 ) one has 

\G (x,y;z)\ < K(x,y;z) = K(x-y,0;z) a.e. x,y, G R 2 

for z sufficiently large negative, K is the Green function of the free Laplacian —A 
in L 2 (R 2 ). Thus h is dominated by ||-0IU 2 (k 2 ) ^ \\ K { Z ) Hl 2 (m 2 ) Wfh^m)- Together 
with the observation Go(x,y; z) = G (y,x; z) (which is a consequence of G (z) being 
a Carleman kernel, see |BGPH theorem 16] and |GMG[ theorem 1]) it implies the 
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claim for large negative z. For any other z G q{Hq), it follows from the first resolvent 
formula 

Rdx,m{z) = Rdx,m( z o) + [Z — Zq) Rq(z) Rdx ,m{zo) . 

□ 

Lemma 3. Consider mappings I m and R mt dx{z) defined as above and let z Q < be 
such that 

(ip,(p) Zo := h(ip,<p) - (zoip,<p) 

is an inner product in D(h). (Recall that h is lower semi-bounded and closed form, 
thus D{h) with (•, -) Zo is a Hilbert space.) Then 1^ = (R m ,d x (zo))* . 

The lemma was proved in [Bj, let us present its proof for the sake of completeness. 

Proof: As both operators I m and Rm,dx{.zo) are bounded, their adjoint operators 
are bounded, too. For any / G L 2 (m) and ip G L 2 (1R 2 ) we have 

W>,-C/)l 2 (R 2 ) = h(R (z )lp,I^f) ~ (Zo-Ro^o)^, C/)x 2 (»2) 
= (Ro(z )tp, I* m f)zo = (ImRo(zo)i>, f)l?{m) 
= (Rm,dx{Zo)lp, f)L 2 (m) = C0> (Rm,dx(z ))* f) L 2^ R 2y 

In the first line, we have used the fact that Ro{zq) is the resolvent corresponding to 
Hamiltonian Ho, associated with h. In the second line, we have simply employed the 
definition of adjoint operator Jj^; since I m maps from D(h) with the inner product 
(•, -) zo to L 2 (m) it reads (0, I* m f)z = {I m <f>, f)i?( m ), € D(h) and / G L 2 {m). □ 

Lemma 4. Assume z G p(Ho) an d f G L 2 (m). Then Rd x .m{z)f G D(h) and 

h{Rdx,m{z)f,ljj) - (zRdx,m(z)f, V , )i 2 (M 2 ) = (fJm^L^m) 

for alii) eD(h). 

Proof: Using the first resolvent formula, one has 

Rdx,m( z )f = Rdx,m( z o)f + [Z — Z )R (z)Rdx !m (z )f, 

where z is the same as in the previous lemma; then the first term equals I^f and 
hence it belongs to D(h). Also the second term belongs to D(h) as the free resolvent 
R {z) maps to D(H ) C D(h). 

To prove the second claim we substitute R dx , m (z)f from the above formula and 
we use lemma EH 

h(R dx>m (z)f,ip) - (zR dx>m (z)f,ip) L 2 ( j S 2 ) = 

= (Rdx,m(z )f,1p) zo - ((Z - Z )Rdx,m(zo)f,4>)L2(R2) 

+ h((z - zo)Ro{z)R d x,m{zo)f, tp) - {z{z - z )Ro{z)R dx ,m(zo)f, ^)l 2 (ir 2 ) 

= (O>V0*> - (0 - Z )R dx , m (z )f, V) L 2 (R 2) 

+ ((Z - Z )R dx , m (z )f,?p)L2(R2) 
= {f, Im^L^im)- □ 
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Theorem 5. Suppose that 1/7 — Rm,m{z) is invertible. Then R(z) given by is 



defined on L 2 (M. 2 ) and it is the resolvent of H, 



Proof: Take arbitrary ip G L 2 (R 2 ) and ip G D(h), then by lemma 0] and Krein's 
formula (jSJ R(z)ip belongs to D{h). We have to check that 

I := h im (R{z)il),(p) - (zRiyz)^^)^^) = (ip,ip) L 2 {R 2y 

Denote g = (1/7 - R^miz))" 1 R m ,d x (z)ip then 

I = h(Ro(z)ip, <p) - (zRo(z)ip, v?)l2(k2) - {I m Ro(z)ip, ^I m ^)^{m) 

+ h(R dXtm (z)g,(p) - (zR dX:m (z)g,ip) L 2 (M 2 ) - {I m Rdx, m {z)g^I m ip) L 2 {m) 

According to lemma |H the fourth and fifth term give together (g, I m (p)L 2 {m), so one 
gets 

/ = 0,y3) L 2 (R2) + ( - R mtdx (z)ip + -g - R mtm {z) g,^ I m ip) , 

V 7 ' L 2 (m) 

finally, employing the definition of g shows that the second term equals zero. □ 
One can fulfil the hypothesis that the operator I/7 — R m ,m{z) is invertibl e easily 
by choosing sufficiently large negative z; it follows from the next lemma, see jBEKS ' 



corollary 2.2]. In the following, ||T|| P)9 denotes the norm of an operator T acting from 
L p (m) to L q (m). 

Lemma 6. There exists z < such that ||7-R m , m (2;)||2 j 2 < 1 for all z < z. 

Proof: Since the measure m belongs to Kato class and 7 is bounded, we can find 
< a < 1 and < b < 00 such that 



\I m ^)(x)\ 2 (l +7(x) 2 )m(cfe) < a h(i/j , i>) + b (i/j, VOi^ 2 ) 

for all ip G D(h). We put z = —b/a, the rhs of the inequality then reads 

a(ip,?p)~ z := ah(ip,4)) - a(zip,ip) L 2 {R 2). 
Next, we take any / G L 2 (m) and z < z and introduce a set 

S z := U; E D(h) : {^) s = l}. 

Consequently, we have 

\I m Rd x , m {z)f{x)\ 2 rf(x)m(dx) < a(R dXim (z)f,R dXtm (z)f) z 
< a sup \(Rd x , m (z)f,ip) z \ 2 = a sup \(f,I m ip)L 2 (m)\ 2 



<a I \f(x) 1 2 m(dx) sup / |/ m ^(x)| 2 m(dx) 

J ip£S z J 



< a 



i2 (m) a(^), = a 2 ||/||i 2(m) □ 



For further use we rewrite the Hamiltonian H im into the form H$ — -//, where 

7m 1 

H —p , Q. —7. . 

J 7m J 7m 

Since both function 7 and measure m are positive, the coupling constant a and 
normalized measure \x are positive, too. The resolvent R(z) = (H jm — z)" 1 then acts 
on arbitrary %p G L 2 (M. 2 ) as follows 

R(z)ip = Ro(z)iP + R dx ,^z)[a - R^z)}" 1 R^ dx (z)^ . (6) 

According to lemma El || -R^^z) 1 1 2,2 is less than 1 for sufficiently large negative z, 
one can prove the same about the norm || ■ || 00,00 using the Kato property of \l. From 
now on we consider only large negative z so that both norms above are less than 1. 
The second term on the rhs of (jUJ) can be substituted by R dXtlM (z)a, where a G L 2 {fi) 
is the unique solution to the equation 

[a - R^(z)]a = R^ dx (z)tp fi - a.e. (7) 

By |BGP1| theorem 16] domain D(H ) is embedded into the space of continuous and 
bounded function on M 2 , thus Ro(z)ip is bounded and continuous in M 2 and the same 
is of course true for function R^d x {z)ip on T. Adding the information about norms 
of we may conclude that a is bounded and continuous on T as well. 



5 Point potentials 

Next consider a magnetic Schrodinger operator Hy a with finitely many point poten- 
tials placed at points a G Y C T, \Y\ denotes the number of potentials. The operator 
is defined via self-adjoint extensions; away from the points from Y it behaves as 
the free operator Hq and the wavefunctions from its domain must have following 
behaviour in the vicinity of each point a EY, 

ip(x) = In \x — a\Lo(ip, a) + Li(ip, a) + 0(\x — a\) (8) 

with coefficients L and L\ fulfilling the boundary condition 

Li(i/;, a) + 27ta(a)L (il),a) = Va G Y. 

For further details concerning point potentials see e.g. [AGHHj for the non-magnetic 
case and |GHSj for the magnetic case. In general, one can choose any real number 
a(a) for each potential independently, here, we make a special choice a(a) = a\Y\ 
for all a EY, with a defined in the previous subsection. 
The resolvent {Hy t a — z)~ X is given by Krein's formula, 

(H Y , a -z)- 1 ^(x)=R (z)^(x)+ [Ha^THy,y')Go(x,y,z)R (z)i;(y'), (9) 

y,y'eY 
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where Ky,a{z) is a matrix \Y\ x \Y\, 

A 



yy 



Y\a-£{y;z) y = y' 
-G (y,y';z) y^y', 



C,(y; z) is the regularized Green function (P). For a homogeneous magnetic field B, we 
can calculate £(a; 2) explicitly; it does not depend on the potential position y and it 



equals -j- 



\B\-z 



2\B\ 



+ 2C E + In ( J|l 



with denoting the Euler constant. The 

second term on the rhs of © can be rewritten as Xlyey Go(x, y; z)q y , where q is the 
\Y I -dimensional vector which solves 



Ro(z)i/>(y) = J2( A rAz))(y,y')%' v ^ Y - 
y'& 



(10) 



One possible way to make matrix A invertible is to take a sufficiently large set Y . 

00, and 



Lemma 7. Let Y n be a sequence of subsets of T with \Y n 
assume that for some a < a we have 



00 as n 



sup 



neN \Y n \ xeY n 



sup 



^2 \G (x,y; z)\ < a < a. 



y£Yn\{x} 

Then there exist C > and ji gN such that A Yn , a ( z ) is invertible and 

1 



y 



< c 



2.2 



for all n> n . 



Here and in the following section, the symbol ||T|| Pjg denotes the operator norm of a 
matrix T, acting from £ P (Y) to £ q (Y) and correspondingly, || ■ \\ p is the norm in £ P (Y). 

Proof: Let us split matrix AY n , a (z) / \Y n \ into the diagonal and non-diagonal part. 
On the diagonal, £(•; z) stays bounded because it is a continuous function on a com- 
pact set T. Therefore the diagonal part behaves as (a + 0{\Y n \~ l ))8 yy i for large n, 
its norm in £ 2 (Y n ) converges to a and so it is invertible. The non-diagonal part, 



denoted by R n , is given by —^G Q (y,y',z)(l 
(see [AGHH , appendix C]) we arrive at 



Syyi). Using Schur-Holmgren bound 



\R 



n 2,2 



< 



\Y 



sup 



n\ x& n 



\G (x,y; z)\ < a < a, 



y£Y n \{x} 



hence the whole matrix is invertible. 



□ 
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6 Approximation 



Now we are ready to formulate the main approximation claim. Its proof follows 
closely the one in |BFT| . 

Theorem 8. Let T be a compact and non-empty set in R 2 and let m be a finite 
positive measure with V := supp(m), which belongs to the Kato class. Let 7 be a 
bounded and continuous function on T, which attains only positive values. Consider 
sufficiently large negative z such that the equation ffi) has a unique solution a with 
a bounded and continuous version on T. Suppose further that there is a sequence of 
sets Y n C T with \Y n \ —>■ 00 as n —>■ 00 and satisfying following three conditions: 

T7F7 E fto -> / f(vM*v) ( n ) 

for any bounded and continuous function f on Y , 

sup r—r sup V] \G (x,y;z)\ < a (12) 

n€N K« x&Y n . , , 

yeY n \{x} 

for some a < a, and finally 

7777 E a (y) G o( x >V> z ) - (i2dx, M (2)<r)(x) 



sup 



yey„\{x} 



(13) 



for n — > 00. T/ien operators Hy n , a converge to H ym in the strong resolvent sense as 
n — > 00. 

Proof: For self-adjoint operators the weak resolvent convergence implies the strong 
resolvent convergence, thus it is enough to prove that 

I n : = (0, (H Yn , a - z)- l ip - {H im - z)~ 1 ip)i^Q i 2) ^ as n -> 00. 

Employing the alternative expressions for both resolvents one gets 



V J/GY„ 



y; z) <? y - RdxA z ) 
where q and a are given as solutions to equations (fTUj) and respectively 

J n = E ^o(^)0(?/) - (R»,dx(z)(t>, 0")i 2 ( M ) 

+ £ i^R0(y) <t(j/) - y L~Mz^{y') a{y') fx(dy'). 



y&. 
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Since Ro(z)<j) has a bounded and continuous version (and it can be identified with 
I^Rq(z)(J)), by hypothesis (fTTj) the last two term vanish in the limit and for the first 
term we only have to show that 



,(n) 



0. 



,(») 



\Y n \q y - °(y)- 

Comparing equations ()10|) and (|7jl. one obtains following expression for i>( n ), 



I»Go(y,y'] z)a{y') n{dy'). 



Last two terms on the rhs vanish because of the hypothesis (|13|) and also the first 
term goes to zero as n —>■ oo since the numerator is a bounded function of y. So when 
we denote the vector with elements given by the rhs as w^ n \ then the norm ||iy^||oo 
tends to zero. 

By hypothesis (|T2*|) and lemma[7[ matrix A.Y n , a {z)/\Y n \ is invertible and the oper- 
ator norm of its inverse in £ 2 (Y n ) is bounded by some C. Hence we can write 



\Y 



since 



< 



1 



\Y„ 



1 



\Y 

- 1 n 



k (n) IU < 



i 



\Y n 



\Y n \ C 



k (n) ||oo 



oo.l 



oo.l 



< \Y n \ || • || 2,2 f° r operators on C' r ™ 



□ 

We have formulated the approximation result for the two-dimensional situation, 
however, it could be proved also in dimension three, provided several modifications 
are made. First of all, the generalized Kato class is different in M 3 . Measure m 
belongs to it if 

f 1 

lim sup / 1 m(dy) = 0, 

£ ^°xGR3 JB(x,e) \X-y\ 

where B(x, e) is the sphere of radius e centred at x. A potential with a zero measure 
support fulfils the condition above only if the codimenion of its support is equal to 
1. Therefore one may approximate for example the potentials supported by compact 
surfaces. 

The main technical difficulty in dimension three comes from the fact that the 
Green function of a magnetic Schrodinger operator has a divergent term which de- 
pends on the given vector potential, see e.g. |BGP2t section 5]. Correspondingly, the 
vector potential also enters the definition of point potentials, hence the definition has 
to be modified as in |KN1| section 4]. Moreover, the divergent term is not logarithmic, 
hence the term In \x — a\ in (jHJ) must be replaced by l/\x — a\. 
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7 Application 



One motivation for this paper was to obtain an alternative method to calculate dis- 
crete spectra. We have proved that any operator H jm , defined as H perturbed by 
measure m with Kato property, can be approximated by point potential Hamiltoni- 
ans in the strong resolvent sense. Hence for each eigenvalue of H im there exists a 
sequence of eigenvalues of the latter operators converging to it. 

It is very natural to apply the approximation to a system whose spectral problem 
is exactly solvable so that we can compare the exact and approximate eigenvalues, 
obtained by numerical calculation. As an essential prerequisite for numerical cal- 
culations is to have an explicit formula for Green function Gq(z), we restrict the 
application only to the case of a homogeneous magnetic field B. Then one can em- 
ploy the expression (J2J). 

To demonstrate how one can use the approximation to calculate discrete spectra, 
let us present following example. Suppose that the potential is attractive and constant 
and it is supported by a circle with radius R. The potential can be thus described by 
two parameters R > and 7 > 0. The easiest choice of point potential operators is the 
following: we place iV points equidistantly along the circle and put a = N/ (2irR'~f) . 
The spectrum of H 1Tn consists of Landau levels 

o c& s{H im ) ={|B|(2m+l): m = 0,l,...} 

and eigenvalues which have split off from the Landau levels because of the presence 
of the potential. Since we consider the potential being attractive, the eigenvalues 
are below the level they have arisen from (and moving further down as the coupling 
constant 7 grows.) 

To find the eigenvalues explicitly, one has to decompose the operator into angular 
momentum subspaces and then to calculate the eigenvalues numerically by solving 
an implicit equation in each subspace, see |ETj . The resulting picture is that there 
is one sequence of eigenvalues in each gap between two adjacent Landau level and 
below the lowest one, with the limit point at the upper Landau level. 

Figure ^ depicts the comparison of the approximate and exact eigenvalues in two 
lowest gaps for two situations which differ only in the coupling constant 7. Figure 
|l(b)| corresponds to a stronger attractive potential 7 = 3, therefore the eigenvalues 



are further from the Landau levels than those in 1(a) where 7 = 1. We observe that 
the approximate eigenvalues tend to the exact ones as the number of point potentials 
grows and that the convergence is slower when the coupling is stronger. One can 
roughly estimate that the convergence rate is of the type cN~ a , where according 
to numerical calculations, a appears to be around 1/2, while coefficient c depends 
strongly on the coupling constant 7. 

A close inspection of the eigenfunctions of the point potential operators would 
reveal that they have logarithmic peaks at the potential sites, since they are given as 
linear combinations of free Green functions, see e.g. |AGHH[ chapter II]. Although 
by |Po21 theorem 3.4], these wavefunctions yield an approximation of wavefunctions 
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Figure 1: The dependence of the approximate eigenvalues on the number of point 
potentials for B = 1, R = 2 and 7 = 1 (a) and 7 = 3 (b). The dashed lines represent 
the exact eigenvalues of H im . 



of H im , the peaks are of course absent in exact eigenfunctions. We believe that this 
fact is partly responsible for the slow convergence of the approximate energies of the 
bound states to the exact ones. 

All the features we have described were observed in the non-magnetic case, see 
[EN2J, with the exception that there one deals only with one gap (— 00, 0) and the 
number of eigenvalues in the gap is finite. In the absence of magnetic field the ground 
state of H im corresponds to angular momentum 1 = 0, the remaining bound states 
correspond to ±Z and they are double degenerate. In magnetic field there is no such 
degeneracy; eigenvalues for angular momenta with opposite signs are different, be- 
cause there is an extra angular momentum coming from the magnetic field. As figure 
|l(b)| suggests, also the approximate eigenvalues (and in particular, their dependence 
on the number of point potentials) behave differently: the eigenvalue crossing the 
Landau level B = 1 tends to the eigenvalue of H im for 1 = 1, while the eigenvalue for 
/ = — 1 is the limit point of the second lowest approximate eigenvalue. 
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